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We present a study of the inclusive production of the X(4140) with the decay to the J/ifcj) final 
state in hadronic collisions. Based on 10.4 fb“^ of pp collision data collected by the DO experiment 
at the Fermilab Tevatron collider, we report the first evidence for the prompt production of df (4140) 
and find the fraction of X(4140) events originating from b hadrons to be fb = 0.39 ± 0.07 (stat) ± 
0.10 (syst). The ratio of the non-prompt X(4140) production rate to the 5° yield in the same channel 
is R = 0.19±0.05 (stat) ±0.07 (syst). The values of the mass M = 4152.5± 1.7(stat)l5’4 (syst) MeV 
and width T = 16.3 ± 5.6 (stat) ± 11.4 (syst) MeV are consistent with previous measurements. 

PACS numbers: 14.40.Cx,13.25.Cv,12.39.Mk 


The X(4140) state [I| was first seen in 2009 as a 
narrow structure in the J/'4’4> system near threshold. 
The CDF Collaboration reported the first evidence [1] 
for this state (then designated y(4140)) in the decay 
5+ -ri X{A140)K~^ -ri J/ipcfK^^ (charge conjugation is 
implied throughout this paper) and measured the invari¬ 
ant mass M = 4143.0 ± 2.9(stat) ± 1.2 (syst) MeV and 
width r = 11.71 |;q (stat) ± 3.7 (syst) MeV. The LHCb 
Collaboration found no evidence for the Ar(4140) state Q 
in a 2.4 standard deviation disagreement with the CDF 
measurement. However, the presence of 1^4140) in 
decay was later confirmed by the CMS Q and DO @ 
Collaborations. The BaBar Collaboration searched for 
resonant production in the J/ipfj) mass spectrum in B'^'^ 
decays and obtained a significance below 2a, but noted 
that the hypothesis that the events are distributed uni¬ 
formly on the Dalitz plot gives a poorer description of the 
data 0. The quantum numbers of the X(4140) state 
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have not been measured. Since both the J/if and (f 
mesons have = 0“1 , the state has positive 

G and C parities. 

A meson decaying into a charmed quark pair might 
be an excited charmonium state. However, the standard 
nonrelativistic quark model of a single cc pair does not 
predict a hadronic state at this mass. Also, at masses 
above the open-charm threshold of 3740 MeV such states 
are expected to decay predominantly to pairs of charmed 
mesons and to have a much larger width than is experi¬ 
mentally observed. It has been suggested that A(4140) 
could be a molecular structure made of two charmed 
mesons, e.g. {Ds,Ds). Other possible states are hybrids 
composed of two quarks and a valence gluon (qqg), four- 
quark combinations (ccss), or states with higher Fock 
components Q . For details see the reviews in Ref. Q and 
Q and references therein. The Belle Collaboration found 
no evidence for A(4140) in the process 77 J/4’4’ 0 , 

making its interpretation as a hadronic molecule with 
spin-parity = 0 "^ or 2 + unlikely. 

In addition to A(4140), the CDF Collaboration re¬ 
ported seeing a second enhancement in the same channel, 
located near 4280 MeV. A similar structure is seen by 
the CMS Collaboration Q at a slightly higher mass of 
4316.7 ± 3.0 (stat) ± 7.3 (syst) MeV. Belle also reports a 
new structure at M = 435 O. 6 I 5 ® (stat) ±0.7 (syst) MeV. 

In this Article we present results of a search for the 
A (4140) resonance in the J/ip4> system produced inclu¬ 
sively in pp collisions, either promptly, by pure QCD, 
or through weak decays of b hadrons. The measured 
production rates are normalized to the rate of the pro¬ 
cess Hg ^ J/4’4’ measured with the same dataset. The 
data sample corresponds to an integrated luminosity of 
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10.4 fb ^ collected with the DO detector in pp collisions 
at 1.96 TeV at the Fermilab Tevatron collider. 

The DO detector consists of a central tracking system, 
calorimeters, and muon detectors [ll|. The central track¬ 
ing system comprises a silicon microstrip tracker (SMT) 
and a central fiber tracker (CFT), both located inside a 
1.9 T superconducting solenoidal magnet. The tracking 
system is designed to optimize tracking and vertexing for 
pseudorapidities | 7 y| < 3, where p = — ln[tan(0/2)], and 0 
is the polar angle with respect to the proton beam direc¬ 
tion. The SMT can reconstruct the pp interaction ver¬ 
tex (primary vertex) for interactions with at least three 
tracks with a precision of 0.004 cm in the plane trans¬ 
verse to the beam direction. The muon detector, posi¬ 
tioned outside the calorimeter, consists of a central muon 
system covering the pseudorapidity region Iryj < 1 and a 
forward muon system covering the pseudorapidity region 
1 < \p\ < 2. Both central and forward systems con¬ 
sist of a layer of drift tubes and scintillators inside 1.8 T 
iron toroidal magnets with two similar layers outside the 
toroids (ll| . 

Events used in this analysis are collected with both 
single-muon and dimuon triggers. Muon triggers require 
a coincidence of signals in trigger elements inside and out¬ 
side the toroidal magnets. Dimuon triggers in the central 
rapidity region require at least one muon to penetrate the 
toroid. In the forward region, both muons are required 
to penetrate the toroid. 

We study a wide range of the J/ip'P invariant mass, 
from threshold to 5.7 GeV, covering both the X(4140) 
and the decay J/ip(j). Candidate events are re¬ 

quired to include a pair of oppositely charged muons in 
the invariant mass range 2.9 < M(^+/i“) < 3.3 GeV, 
consistent with J/V' decay, accompanied by two ad¬ 
ditional particles of opposite charge, assumed to be 
kaons, with pT > 0.4 GeV and 1.011 < M{K'^K~) < 
1.030 GeV. In the event selection, both muons are re¬ 
quired to be detected in the muon chambers inside the 
toroidal magnet, and at least one of the muons is re¬ 
quired to be also detected outside the iron toroid [l^ . 
Each muon candidate is required to match a track found 
in the central tracking system, and each of the four final- 
state tracks is required to have at least one SMT hit 
and at least one GET hit. The dimuon invariant mass is 
constrained to the world-average J/4’ mass [I|, and the 
four-track system is constrained to a common vertex. To 
reconstruct the primary vertex, tracks are selected that 
do not originate from the candidate, and a con¬ 

straint is applied to the average beam position in the 
transverse plane. We require candidates to have 

5 < pt < 20 GeV and rapidity \y\ < 2. 

We define the signed decay length of the J/4’4> system, 
Lxy, to be the vector pointing from the primary vertex 
to the decay vertex, projected onto the direction of the 
transverse momentum. The distribution of L^y for the 
selected events is shown in Eig.[T] 

We focus on two ranges of the invariant mass, 

< 4.36 GeV and 4.8 < < 5.7 GeV. 



FIG. 1: Transverse decay-length distribution of J/iKj) candi¬ 
dates. The vertical lines define the three regions discussed in 
the text. 


The low-mass range includes the V(4140) state. The 
high-mass range includes the reference decay process, 
Bg J/ipip. Background arises primarily from non¬ 
resonant pairs in the (p mass window. At low Lj.y, back¬ 
ground comes from J/ip mesons directly produced in pp 
collisions combined with random particles from the un¬ 
derlying event. At higher values of background con¬ 
sists of J/ip mesons paired with random products of b 
hadron decays. 

We divide the data in each mass range into three in¬ 
dependent subsamples according to the value of L^y'. (1) 
—0.025 < Lxy < 0 cm, (2) 0 < L^y < 0.025 cm, and (3) 
Lxy > 0.025 cm. Region 1 includes half of the prompt 
events and almost no R-decay events (the fit result shown 
in Table U is 37 ± 26 events). Region 2 includes the re¬ 
maining half of all prompt events and a fraction of non¬ 
prompt events. The rest of the non-prompt events pop¬ 
ulate region (3). Given the average resolution of 0.006 
cm in L^y^ we assume that the fraction of prompt events 
in Region 3 is negligible. We perform binned maximum 
likelihood fits to the distributions of the invariant 

mass for events in the six subsamples defined above. In 
the fits in the mass region, the signal is described by 
a Gaussian function and background is described by a 
second-order Chebychev polynomial. We also allow for 
the presence of the decay where we set the 

mass to the world-average B^ mass, and we find no evi¬ 
dence of a signal. The fit for Region 3 yields 3166 ± 81 
events. 

In fitting the low mass range, we assume a signal de¬ 
scribed by an tS-wave relativistic Breit-Wigner function 
convolved with a Gaussian resolution of a{M) = 4 MeV. 
The background is parametrized by the function f{m) oc 
m ■ {rn ^where mthr is the kinematic 
threshold, and ci and C 2 are free parameters. Eor events 
in the L^y Region 3, we allow the signal mass and width 
parameters to vary. The fit yields 616± 170 signal events, 
a mass of 4152.5 ± 1.7 MeV, and a width of 16.3 ± 5.6 
MeV. The statistical significance of the signal, based on 
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the increase of the likelihood with respect to the fit with 
no signal, —2Aln>C = 42.5 for 3 degrees of freedom, is 
5.9 standard deviations. For the fits in L^y Regions 1 
and 2 we set the mass and width to the Region 3 values. 

The mass distributions with superimposed fits for both 
mass regions and for all three L^y are shown in Fig. [2] 
The A(4140) and 5° yields are presented in Table HI We 
also show the expected number of A(4140) events origi¬ 
nating from 6 -hadron decays in the two low L^y regions 
assuming that the Lxy distribution of the “non-prompt” 
X(4140) is similar to that of 5°. For the Regions 1 and 
2 , we find an excess of signal events, indicating prompt 
production of A(4140). For events in Region 2, the in¬ 
crease in the likelihood between the fit with a free signal 
yield and the fit with the expected non-prompt contribu¬ 
tion only, —2Aln£ = 23.6, corresponds to a statistical 
significance of 4.9cr for the net prompt signal. The sta¬ 
tistical significance of the total signal in this Lxy region 
is 6 . 2 ( 7 . For Region 1, the corresponding values of sta¬ 
tistical significance are 3 . 9(7 and 4 . 2 ( 7 . If the mass and 
width parameters are allowed to vary, the fit for Region 
2 gives the total yield N = 932 ± 216, M = 4146.8 ± 2.4 
MeV, and F = 15.8 ± 3.8 MeV. The data in Region 1 do 
not yield a stable fit. Fixing the A(4140) mass to 4152.5 
MeV in this region, as obtained in Region 3, we fit a total 
yield oiN = 601 ± 205 and F of 19.8 ± 5.9 MeV. 

There are several uncertainties that may affect mea¬ 
surements of the X(4140) yield, mass, and width, the 
ratio R of the yields of non-prompt V(4140) and 
and the fraction fy of all V(4140) events that originate 
from weak decays of b hadrons. 

The mass resolution of 4.0 MeV, obtained in simula¬ 
tions, is in agreement with an approximately linear rise 
with the released kinetic energy for decays with a similar 
topology: tjj{2S) —>■ J/, V(3872) —>• J/Vctt+tt”, 
and the decay 5° —>■ J/ipcj). We assign an uncertainty of 
±0.1 MeV to the resolution at the V(4140) mass. 

We assign an asymmetric uncertainty of 3 MeV to the 
J/ipfj) mass scale in the vicinity of V(4140), based on the 
range of the mass deficit between 1 MeV and 5 MeV com¬ 
pared to world-average values, found in several channels 
with this topology. We assign the uncertainty in the sig¬ 
nal model, taken from the range of results obtained with 
relativistic and nonrelativistic Breit-Wigner shapes and 
a relativistic P-wave Breit-Wigner shape. Simulations 
show the event reconstruction and selection efficiency to 
be independent of the M(J/±^) invariant mass, with a 
possible variation of ±10% Q. The possible variation 
of the efficiency within the V(4140) mass range affects 
the mass, width, and yield of the signal. To assess the 
effects of the fitting procedure and background size and 
shape, we vary the fitting mass range and bin size. Some 
of the single-muon triggers include a trigger term requir¬ 
ing a presence of tracks with non-zero impact parameter. 
Events recorded solely by such triggers constitute approx¬ 
imately 5% of all events. Assuming that such triggers are 
100 % efficient for events originating from weak decays of 
b hadrons and reject all prompt events, we apply a 5% 


correction to the prompt yield. We assign a systematic 
uncertainty of ±5% on the fraction /{, due to this correc¬ 
tion. Finally, our assumption of the equality of the rela¬ 
tive rates in regions (1) - (3) for the non-prompt A(4140) 
and P° is based on expectation of the equality of the av¬ 
erage lifetime of 6-hadron parents of the A(4140) and 
that of the in the J/ip(j) channel. The world-average 
of the Pg lifetime is 6% lower than the lifetime averaged 
over all 6 hadron species [l|. We assign an asymmetric 
uncertainty in the ratio R and the fraction fy based on 
this difference. The systematic uncertainties are summa¬ 
rized in Table HU 

We test the stability of the results to the event selection 
by changing the (j> mass window to 1.012 < M{K'^K~) < 
1.029 GeV. As additional cross-checks, we perform fits to 
subsamples corresponding to the transverse momentum 
ranges 5 < pr < 10 GeV and 10 < pt < 20 GeV; to early 
and late data-taking periods; and to events in the central 
(IpI < 1) and forward rapidity regions. In each case, the 
background shape in the two subsamples is well described 
by the same functional form although it requires different 
values of the parameters. In all cases the sums of the 
resulting signal yields agree with the total yield within a 
few events. 

Our measured values of the mass and width of the 
X (4140) state are compared with earlier measurements in 
Table Hm The ratio of the A(4140) to B^ yield for events 
with Lxy > 0.025 cm is R = 0.19±0.05 (stat)±0.07 (syst). 
After correcting for the efficiency of this Lxy cut and 
for the trigger bias, we find the fraction of A(4140) 
events originating from 6 hadrons to be ft = 0.39 ± 
0.07 (stat) ± 0.10 (syst). The yield for the A(4140) state 
at Lxy > 0.025 cm can also be compared with the yield 
of 52 ± 19 events of A(4140) from the decay process 
P+ —>■ J/ipcpK'^ obtained by DO for the same data 
set. After correcting for a factor of 2.5 ± 0.5 for the ef¬ 
ficiency of the full reconstruction of the B^ decay and 
lower kaon px threshold, we expect the yield from the 
B'^ decay to be ~ 130 ± 60 events in this analysis. Our 
observed yield of 616 ± 170 events exceeds this estimate 
suggesting that decays of 6 hadrons other than P+ con¬ 
tribute to the non-prompt production of A(4140). 

The invariant mass distributions presented in 

Fig. [2] show no evidence for states in the mass region 
4250 < M{J/'ijj(j)) < 4375 MeV. Fits allowing for the 
the states reported by CDF, CMS, and Belle, at Lxy > 
-0.025 cm, yield 267 ± 276, -283 ± 468, and -325 ± 254 
events, respectively. Using the CLg method [l^, we ob¬ 
tain the 95% upper limits of 744, 557, and 338 events 
for the three states. In this upper limit calculation we 
did not account for systematic uncertainties as they were 
checked to have a negligible impact. The corresponding 
95% upper limits on the rates relative to the total yield 
of the A(4140) state are 0.43, 0.31, and 0.19. 

In summary, we have carried out the first search for in¬ 
clusive production of the state A(4140) in hadronic col¬ 
lisions. We find strong evidence for its direct, prompt 
production, and observe its production in weak decays of 
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FIG. 2: (color online) Invariant mass distribution of candidates in the mass window around (left) and (right) X(4140), 
for events with (a,b) —0.025 < L™ < 0 cm, (Cj^ 0 < L^y < 0.025 cm and (e,f) L^y > 0.025 cm. The arrows indicate the 
structures seen by CDF Q], CMS and Belle [1^. The signal and background models are described in the text. 


TABLE I: Summary of event yields in three L^y regions and their sum for and X(4140). For Regions 1 and 2 the mass 
of X(4140) is assumed to be 4152.5 MeV and the width is taken to be 16.3 MeV. Also shown are the deduced yields for the 
non-prompt and prompt production of -A(4140). The uncertainties are statistical. 


Parent 

— 0.025 < Lxy < 0 cm 

0 < Lxy < 0.025 cm 

Lxy > 0.025 cm 

Sum 

B'i 

191 ± 143 

804 ± 169 

3166 ±81 

4161 ± 236 

X(4140) 

511 ± 120 

837 ± 135 

616 ± 170 

1964 ± 248 

A(4140) non-prompt 

37 ±26 

156 ± 54 

616 ± 170 

809 ± 175 

A(4140) prompt 

474 ± 123 

681 ± 149 

= 0 

1155 ± 193 
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TABLE II: Summary of systematic uncertainties. 


Source Mass (MeV) Width (MeV) Rate non-prompt (%) Rate prompt (%) 

Mass resolution 

±0.1 

±0.2 

±1 

±1 

Mass bias 

+3 

-0 

- 

- 

- 

Efficiency 

±4 

±5 

±4 

±4 

Signal model 

±1 

±2.7 

±13 

±15 

Fitting range 

±3 

±7.0 

±20 

±6 

Bin size 

±1.6 

±7.0 

±25 

±10 

Trigger bias 

- 

- 

- 

±5 

Mean lifetime 

- 

- 

-1.5 

±1.5 

Total 

- 5.4 

±11.4 

±35 

±19 


TABLE III: Summary of A(4140) measurements. 


Experiment 

Process 

Mass (MeV) 

Width (MeV) 

CDF 

CMS ]£ 

B+ J/i)(j)K+ 

B+ J/ip(l>K+ 

4143.0 ±2.9 ±1.2 
4148.0 ±2.4 ±6.3 

11.7+50 ±3.7 

28t)5 ± 19 

DO [a 

B+ J/il}4>K+ 

4159.0 ±4.3 ±6.6 

19.9 ± 12.618.0 

DO (this work) pp —> J/'ip4> ± anything 

4152.5 ±1.7+5;4 

16.3 ±5.6 ±11.4 


b hadrons with a rate exceeding the expected rate for the 
known decay —>■ J/'ip(j)K'^. The significance of the 

prompt production, including systematic uncertainties, 
is 4.7cr. This is the first evidence for the prompt produc¬ 
tion of AT (4140). The significance of the non-prompt pro¬ 
duction, including systematic uncertainties, is 5.6 ct. The 
non-prompt production rate of X(4140) relative to 
observed in the same final state is i? = 0.19±0.05 (stat)± 
0.07 (syst). Assuming a relativistic Breit-Wigner line 
shape, we measure the mass and width of the X(4140) 
state to be M = 4152.5 ± 1.7(stat)l54 (syst) MeV and 
width r = 16.3 ± 5.6 (stat) ± 11.4 (wst) MeV, consistent 
with previous measurements 
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